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ABSTRACT 


There  has  been  recent  interest  in  the  use  of  quadtrees 
to  represent  regions  in  an  image.  It  thus  becomes  desirable 
to  levelop  efficient  methods  of  conversion  between  quad¬ 
trees  and  other  types  of  region  reproar-ntations.  This 
paper  presents  an  algorithm  for  converting  from  quadtrees 
to  a  simple  class  of  boundary  codes. 
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Circular  Model 
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3.  1 

Wo  will  first  compute  IQ  and  fj  for  circles. 
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Considering  a  chord  passing  through  a  point  at  a  distance 
r  from  the  center  of  a  circle  of  radius  R,  r  =-R  (Kig.  5),  we  hav 
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This  type  of  boundary  representation  is  railed  a  chain  code. 
>neralized  chain  codes,  involving  more  than  four  directions, 
can  also  be  used.  Chain  codes  provide  a  very  compact  region 
representation,  and  make  it  easy  to  detect  features  of  the  reqion 
boundary,  such  as  sharp  turns  (“corners")  or  concavities.  On 
the  other  hand,  it  is  harder  to  determine  properties  such  as 
elonqatedness  from  a  chain  code,  and  it  is  also  difficult  to  per- 
form  operations  such  as  union  and  intersection  on  regions  repre¬ 
sented  by  chain  codes.  A  general  introduction  to  chain  codes  and 
their  uses  can  be  found  in  (11. 

Another  class  of  region  representations  involves  various  types 
of  maximal  “blocks*  that  are  contained  in  a  given  reqion.  For 
example,  we  can  represent  a  region  R  as  a  linked  list  of  the  runs 
(of  pixels)  in  which  R  meets  the  successive  rows  of  the  array  12) . 
Here  each  “block"  is  a  1-by-m  rectangle,  where-  m  is  the  run 
length;  the  runs  are  the  largest  such  blocks  that  R  contains, 
and  R  is  determined  by  specifying  the  initial  j>ointR  (or  centers) 
and  lengths  of  the  runs.  Alternatively,  we  can  represent  R  by 
the  set  of  maximal  square  blocks  (or  blocks  of  any  other  desired 
shape)  that  it  contains;  here  R  is  determined  by  specifying  the 
centers  and  radii  of  these  blocks.  This  representation  is  called 
the  medial  axis  transformation,  or  MAT  (3).  It  is  somewhat  less 
compact  than  chain  code  |4),  but  it  has  advantages  with  respect 
to  performing  union  and  intersection  operations  or  detecting 


properties  such  as  elongatedness  (in  terms  of  the  smallness  of 
the  radii  relative  to  the  number  of  centers). 

There  has  been  recent  interest  in  an  approach  to  region  rep¬ 
resentation  based  on  successive  subdivision  of  the  array  into 
quadrants.  If  the  region  does  not  cover  the  entire  array,  we  sub¬ 
divide  the  array,  and  repeat  this  process  for  each  quadrant,  each 
subquadrant,...  as  long  as  necessary,  until  we  obtain  blocks 
(possibly  single  pixels)  that  are  entirely  contained  in  the 
region  or  entirely  disjoint  from  it.  The  resulting  blocks  for 
the  region  of  Figure  la  are  shown  in  Figure  lb.  This  process  can 
be  repres  nted  by  a  tree  of  degree  4  (for  brevity:  a  quadtree) 
in  which  the  entire  array  is  the  root  node,  the  four  sons  of  a 
node  are  its  quadrants,  and  the  leaf  nodes  correspond  to  those 
blocks  for  which  no  further  subdivision  is  necessary.*  The  quad¬ 
tree  representation  for  Figure  lb  is  shown  in  Figure  lc.  Note 
tnat  here  again  we  are  representing  the  region  as  a  union  of 
maximal  blocks,  but  this  time  the  blocks  must  have  standard  sizes 
and  positions  (powers  of  2).  Since  the  array  was  assumed  to  be 
2n-by-2n,  the  tree  height  is  at  most  n.  This  method  of  region 
representation  was  proposed  by  Klinger  16-7);  it  has  also  been 
used  for  image  representation  (e.g.,  (8-11)).  It  is  relatively 

compact,  and  is  also  well  suited  to  operations  such  as  union  and 
intersection,  and  to  detecting  various  reuion  properties.  A 


•  The  quadtree  region  representation  described  here  should  not 
be  confused  with  the  quadtree  representation  of  two-dimensional 
point  data  introduced  by  Finkel  and  Bentley  ( S ] . 


recent  Ph.D.  thesis  by  Hunter  [12]  in  the  domain  of  computer 
graphics  develops  a  variety  of  algorithms  for  the  manipulation 
t  ;uadtiee  region  representations.  Those  algorithms,  however, 
a  1  low  a  node  to  store  the  list  of  coordinate  points  that  describe 
the  polygon  trom  which  the  quadtree  was  constructed. 

Since  the  quadtree  and  border  representations  both  have  eompu 
tational  advantages,  it  is  of  interest  to  develop  methods  of 
converting  from  one  representation  to  the  other.  We  shall  now 
present  an  algorithm  for  deriving  a  clockwise  boundary  code  from 
the  juadtree  ref  resentation  of  .1  given  region. 


2. 


Definitions  and  notation 


Let  each  node  in  a  quadtree  be  stored  as  a  record  containing 
six  fields.  The  first,  five  fields  contain  pointers  to  the  node's 
father  and  its  four  sons,  labeled  NW,  NE,  SW,  and  SE.  These 
items  will  be  referenced  as  FATHER (P) ,  NWSON(P),  NESON(P),  SWSON(P), 
and  SESON(P),  respectively,  for  a  given  node  P.  The  sixth  field, 
named  NODETYPE,  describes  the  contents  of  the  block  of  the  array 
which  the  node  represents,  i.e.,  WHITE  if  the  block  contains  no 
pixels  in  the  region,  BLACK  if  the  block  contains  only  pixels  in 
the  region,  and  GRAY  if  it  contains  pixels  of  both  types. 

Let  the  four  sides  of  a  node's  block  be  called  its  N,  S,  E, 
and  W  sides.  Two  nodes  are  said  to  be  adjacent  along  the  N  side 
of  the  first,  for  example,  if  the  pair  of  blocks  represented  by 
these  n‘  des  toucn  along  tnat  side  (not  just  at  a  corner) .  Given  this 
notation  for  a  node's  sides  and  quadrants,  as  illustrated  in 
Figure  2,  we  now  define  several  functions  which  conveniently 
describe  the  geometry  of  these  labelings.  Throughout  this  paper, 
we  use  T  and  U  as  side  variables,  K  and  L  as  quadrant  variables, 
and  P,  Q,  X,  Y,  and  Z  as  node  variables. 

Define  CSIDE(T)  to  be  the  side  which  is  adjacent  to  side  T 
in  the  clockwise  direction,  e.g.,  CSIDE('N')  ■  'E*.  Similarly, 
define  CCSIDE(T)  and  OPSIDE(T)  to  be  the  sides  which  are  adjacent 
in  the  counterclockwise  direction  and  on  the  opposite  side  from 
T,  respectively.  For  example,  CCSIDE('N')  »  *W'  and  OPSIDE('N') 

■  'S'.  The  value  of  the  boolean  function  AD J ( T , K )  is  true  if 


and  only  if  quadrant  K  is  adjacent  to  side  T  of  the  node's  block, 
e.g.,  ADJ ( 'N* , 'NW' )  -  true. 

Define  REFLECT (T, K)  to  be  the  quadrant  in  which  quadrant  K 
lies  after  beiny  reflected  about  the  side  T  axis.  For  example, 
REFLECT ( ' W ' , ' NE ' )  -  'NW'  and  REFLECT (' N ',' NW • )  -  'SW*.  QUAD(T,U) 
is  defined  to  be  the  quadrant  which  touches  the  corner  formed  by 
sides  T  and  U  (if  T  and  U  are  opposite  sides,  the  function  is 
undefined),  e.g.,  QUAD ( ' E' , ' N ' )  -  ' NE ' .  LINK (T)  returns  the 
chain  code  direction  associated  with  side  T.  If  we  are  using 
the  coding  described  in  Section  1,  this  implies  LINK('N')  *  0, 

LINK ( ' W' )  »  1,  LINK ('S')  -  2,  and  LINK ( ' E* )  *  3. 

It  is  also  convenient  to  be  able  to  access  certain  properties 
of  a  node  concisely.  Therefore,  define  SONTYPE(P)  to  be  the 
quadrant  type  of  P  relative  to  P's  father,  e.g.,  SONTYPE(P)  ■  'NW' 
iff  NWSON (FATHER (P) )  -  P.  Let  SON(P,K)  be  the  Kth  son  of  node  P, 
where  K  is  a  quadrant.  For  example,  SON(P,'NK')  ■  NWSON (P). 
Finally,  define  DEPTH (P )  to  be  the  length  of  the  path  from  node 
P  to  the  root  of  the  quadtree,  TREEDEPTH  to  be  the  maximum  depth 
of  any  node  in  the  quadtree,  and  LEVEL(P)  to  be  TREEDEPTH-DEPTH (P) . 


2. 


Definitions  and  notation 


Let  each  node  in  a  quadtree  be  stored  as  a  record  containing 
six  fields.  The  first,  five  fields  contain  pointers  to  the  node's 
father  and  its  four  sons,  labeled  NW,  NE,  SW,  and  SE.  These 
items  will  be  referenced  as  FATHER (P)  ,  NWSON(P),  NESON(P),  SWSON(P), 
and  SESON(P),  respectively,  for  a  given  node  P.  The  sixth  field, 
named  NODETYPE,  describes  the  contents  of  the  block  of  the  array 
which  the  node  represents,  i.e.,  WHITE  if  the  block  contains  no 
pixels  in  the  region,  BLACK  if  the  block  contains  only  pixels  in 
the  region,  and  GRAY  if  it  contains  pixels  of  both  types. 

Let  the  four  sides  of  a  node's  block  be  called  its  N,  S,  E, 
and  W  sides.  Two  nodes  are  said  to  be  adjacent  along  the  N  side 
of  the  first,  for  example,  if  the  pair  of  blocks  represented  by 
these  n> des  toucn  along  tnat  slue  (not  just  at  a  corner).  Given  tnis 
notation  for  a  node's  sides  and  quadrants,  as  illustrated  in 
Figure  2,  we  now  define  several  functions  which  conveniently 
describe  the  geometry  of  these  labelings.  Throughout  this  paper, 
we  use  T  and  U  as  side  variables,  K  and  L  as  quadrant  variables, 
and  P,  Q,  X,  Y,  and  Z  as  nod*  variables. 

Define  CSIDE(T)  to  be  the  side  which  is  adjacent  to  side  T 
in  the  clockwise  direction,  e.g.,  CSIDE('N')  «  'E'.  Similarly, 
define  CCSIDE(T)  and  OPSIDE(T)  to  be  the  sides  which  are  adjacent 
in  the  counterclockwise  direction  and  on  the  opposite  side  from 
T,  respectively.  For  example,  CCSIDE('N')  -  *W'  and  OPSIDE('N') 

■  'S'.  The  value  of  the  boolean  function  ADJ(T,K)  is  true  if 


For  example,  the  three  cases  shown  in  Figure  3  are  treated  as 
fol lows: 


a)  In  this  case  X  cannot  extend  beyond  P,  though  it  may  be 
bigger  than  Q.  If  X  is  black,  the  new  pair  is  (X,Q) ,  the 
boundary  turns  left,  and  the  new  link  is  the  west  side 

of  Q  (if  X  is  larger  than  Q)  or  the  east  side  of  X 
(otherwise).  If  X  is  white,  the  new  pair  is  (P,X),  the 
boundary  does  not  turn,  and  the  new  link  is  the  north 
side  of  X.  See  Figure  4a-b. 

b)  In  this  case  X  cannot  extend  beyond  Q,  thouqh  it  may  be 
bigger  than  P.  If  X  is  black,  the  new  pair  is  (X,Q),  the 
boundary  does  not  turn,  and  the  now  link  is  the  south 
side  oi  X.  If  X  is  white,  the  new  pair  is  (P,X),  the 
boundary  turns  right,  and  the  new  link  is  the  west  side 
of  P  (if  P  is  no  larger  than  X)  or  the  east  side  of  X 
(otherwise).  See  Figure  4c-d. 

c)  Assume  that  the  region  is  4-connected,  so  that  blocks 
touching  only  at  a  corner  are  not  regarded  as  adjacent. 
Then  if  both  X  and  Y  are  black  (they  need  not  be  distinct 
nodes),  the  new  pair  is  (Y,Q) ,  the  boundary  turns  left, 
and  the  new  link  is  the  east  side  of  Y(if  X  and  Y  are 
distinct  and.  Y  is  no  bigger  than  Q)  or  the  west  side  of 

0  (otherwise;  note  that  if  X  and  Y  are  the  same  node, 

Y  must  extend  at  least  to  the  end  of  Q  because  neighboring 
nodes  can't  properly  overlap.)  If  X  is  black  and  Y  is 
white,  the  new  pair  is  (X,Y),  the  boundary  does  not  turn. 


and  the  link  is  the  south  side  of  X  (if  X  is  no  larger 
than  Y)  or  the  north  side  of  Y  (otherwise) .  If  X  is 
white,  the  new  pair  is  (P,X)  regardless  of  the  node 
type  of  Y  by  virtue  of  the  4-connectedness  property. 

The  new  link  is  the  west  side  of  P  (if  P  is  no  bigger 
than  X)  or  the  east  side  of  X  (otherwise);  the  boundary 
turns  right.  See  Figure  4e-i.  Alternatively,  we  could 
specify  a  similar  set  of  riles  if  we  considered  the 
region  to  be  8-connected. 


4 .  Formal  statement  of  the  algorithm 

The  following  ALGOL-like  procedures  specify  the  complete 
boundary  following  algorithm.  The  main  program  finds  and  marks 
the  initial  (black,  white)  pair  of  nodes.  NEXTLINK  is  the 
recursive  algorithm  which  outputs  the  link  associated  with  the 
current  pair,  and  then  finds  the  next  pair.  Algorithm  F1ND_ 
NEIGHBOR  is  used  to  find  the  adjacent  nodes  X  and  Y.  OVERLAP 
determines  which  of  the  three  cases  shown  in  Figure  3  applies 


to  the  current  pair. 


begin 


comment  given  quadtree  QUADTREE,  find  a  black  node  P  which 

is  on  the  region  boundary  and  has  no  black  not'  *s  adjacent 
to  its  north  side.  Q,  a  white  node,  is  the  westerr-mr  st 
of  P’s  northern  neighbors.  The  pair  (P,Q)  defines  hi 
initial  chain  segment; 
node  P , Q ; 

P  •  ROOT  (QUADTREE); 

whl le  NODETYPE (P)  -  ’GRAY’  do 

if  NODET Y  P  E  ( NWSON ( P ) ) t  ’WHITE’ 
then  P- NWSON (P) 

e  l  se  i  f  NODETY PE  ( NESON ( P ) ) f ’ W»  I TE ’ 
then  P- NESON  (P) 

else  if  NODETYPE  (SWSON (P) ) WHITE* 
then  P* SWSON (P) 
else  P* SESON (P) ; 

Q  -  FIND  NEIGHBOR  (P,’N’,’W’)j 
Mark  (P,Q)  as  starting  pair; 

NEXT  LINK  (P , Q, *N’ ) 


end 


procedure  NEXT  LINK  (P,Q,T) : 
begin 

comment  given  adjacent  leaf  nodes  P,  Q,  where  P  is  black,  Q  is 
white,  and  side  T  of  P  touches  0,  output  the  chain 
description  for  this  section  of  the  boundary; 

node  P,Q,X,Y,Z;  side  T;  integer  i; 

if  (P,0)  is  the  marked  starting  pair  and  we've  seen  it 
before  then  halt; 


for  1*1  step  1  until 


,MIN (LEVEL (P) , LEVEL (0) ) 


do  pirint  (LINK  (T)  )  ; 


comment  determine  next  pair  of  nodes  and  their  common  : ;de 
and  recursively  call  NEXT  LINK; 

Z ‘OVERLAP  (P.Q.  CSIDE(T)); 

If  Z-P 

then  begin 

comment^  Black  overlaps  white; 

X* FIND  NEIGHBOR  (0.  CSIDE(T).  OPSIDE(T)); 
if  NODETYPE (X)  -  'WHITE' 

then  NEXT  LINK  <P,X,T)  /‘Figure 

else  NEXt  LINK  <X,Q,  CCSIDE(T))  /‘Figure 


4a‘/ 

4b‘/ 


end 


else  if  Z»Q 

then  begin 

comment  white  overlaps  black; 

X* FIND  NEIGHBOR  (P,  CSIDE(T).T)} 
if  NuDETYPE (X)  -  ’WHITE’ 

then  NEXT  LINK  (P,X,  CSIDF.(T))  /‘Figure  4c*/ 
else  NEXT  LINK  (X.Q.T)  /‘Figure  4d*/ 

end 

else  begin 

comment  Black  and  white  aliqncd; 

X*  FIND  NEIGHBOR  (P,  CSIDE(T).T); 

If  NODETYPE (X)  -  ’WHITE* 

then  NEXT  LINK  (P,X,  CSIDE(T))  /‘Figs.  4e  and  4q*/ 
else  begin 

Y* FIND  NEIGHBOR  (0,  CCS IDE (T) ,T) ; 
if  NODETYPE (Y)  -  ’BLACK’ 

then  NEXT  LINK  (Y,Q,  CCS  IDF.  T)  )  /•Fios.  4',4i*/ 
else  NEXT  LINK  (X.'.T)  /*Fiqure  4h*/ 

end 

end 


end 


procedure  FIND  NEIGHBOR  (P,T,U) : 

begin 

comment  Given  node  P,  return  node  Q  which  is  adjacent  to  side 
T  of  P  and  touches  P's  TU  corner; 
node  P.Q,  side  T,U; 

Q*P; 

STACK* empty 

while  ADJ (T , SONTYPE (0) )  do 
begin 

comment  find  the  nearest  common  ancestor  of  P  and  0; 
push  SONTYPE (Q)  onto  STACK; 

0* FATHER (0) ; 
end; 

0* SON  (FATHER {Q) .  REFLECT  (T,  SONTYPE  (.,)); 
while  NODETYPE (0)  -  'GRAY'  do 
begin 

comment  follow  reflected  path  back  down  tree  to  locate 

neighbor.  If  STACK  is  empty,  then  0  is  smaller  than  P; 
i f  STACK  not  empty 

then  be  jin  0*SON  (0,  REFLECT (T ,  top (STACK) )) ; 
pop  STACK 

end 

else  0* SON  (Q.  QUAD (OPS IDE (T) ,U) ) 

end; 

return  (Q) 


end 


procedure  OVERLAP  (P,Q,T) : 

begin 

comment  given  two  nodes  P  and  Q  which  touch  along  a  side  U,  where 
CSIDE(U)*»T,  determine  whether  P  extends  farther  in  the 
T  direction  than  Q  (return  P) ,  Q  extends  farther  than  P 
(return  Q) ,  or  both  their  Tth  sides  are  aligned 
(return  0) : 

node  P,U,Lo,Hi;  side  T;  integer  Dp,  Dq,  Diff,  i; 

Dp  •  DEPTH (P ); 

Dq  -  DEPTH (Q) ; 


if  Dp“Dq  then  return  (0)  ; 
if  MAX (Dp.Dq)  -  Dq 

then  begin  Lo*P;  Hi*Q  end  else  begin  Lo»Q;  Hi*P  end 
Dlff*  Dp-Dq 

comment  The  smaller  of  tie  two  nodes  cannot  extend  farther  than 
the  other  because  this  would  imply  that  P  and  0 
properly  overlap,  which  is  impossible.  At  best  the 
smaller  node  can  be  aligned  with  the  other  one,  and  th's 
occurs  if  and  only  if  the  smaller  node  is  at  the  extreme 
T  side  relative  to  the  nearest  common  anc-  stor  of  P  and  0 
for  i * 1  step  1  unti 1  D l f  f  do 

begin 

^f  ADJ  (T .SONTYPE (Lo) )  then  return  (Hi); 

Lo*FATHER  (Lo) 


end; 

return  (0) 


end 


5.  An  example 


Let  us  consider  the  quadtree  shown  in  Figure  1.  Notice 
that  this  quadtree  has  57  n«>>i.-  whereas  an  array  representation 
would  have  required  a  16  by  16  (-2‘j6-cell)  logical  array.  The 
main  program  finds  the  initial  pair  (19,13)  and  then  calls 
NEXT  LINK  to  do  the  boundary  following  from  this  point.  Table 
1  specifies  the  arguments  of  NEXT  LINK  at  each  recursive  call 
along  with  the  link  output  as  a  result  of  this  pairing. 
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Sequence  of  calls  of  NEXT  LINK  for  quadtree  in  Figure  1  . 
The  exponent  on  the  output  indicates  the  number  of  unit 
len  :th  links  associated  with  the  boundary  specified  by 
the  jiven  (P,Q)  [air. 


Table  1  . 


6 .  Analysis 


The  speed  of  this  boundary  coding  algorithm  is  determined  by 
the  procedure  NEXT  LINK  which  is  called  for  each  boundary  segment 
associated  with  a  (black,  white)  adjacent  node  pair.  The  time 
required  to  output  the  individual  links  is  proportional  to  the 
region’s  perimeter,  where  perimeter  is  defined  to  be  the  number 
of  unit-square  pixels  on  the  region's  boundary  (not  the  number 
of  black  nodes  which  are  adjacent  to  white  nodes).  The  additional 
cost  of  the  algorithm  is  determined  by  the  calls  to  procedures 
OVERLAP  and  FIND  NEIGHBOR.  Given  a  pair  (P,Q)  of  adjacent  black, 
white  nodes,  each  procedure  is  called  once  in  order  to  find  the 
next  pair.  The  time  requirt  J  1  measured  by  the  number  of  nodes 
visited  and  depends  on  the  relative  positions  of  P  and  j,  and 
adjacent  nodes  X  and  Y  in  the  quadtree. 

Procedure  OVERLAP  takes  time  proportional  to  the  sum  of  the 
depths  of  nodes  P  and  Q  in  the  quadtree,  althouqh  this  can  be 
reduced  to  the  sum  of  the  path  lengths  to  their  nearest  common 
ancestor.  Further  refinements  can  be  made  based  on  the  relation¬ 
ship  between  P  and  0  so  that  OVERLAP  doesn't  need  to  be  called 
for  each  call  to  NEXT  LINK.  For  example,  having  determined  that 
P  extends  beyond  Q  and  X  is  smaller  than  Q  (which  we  can  detect 
while  searching  for  j()  ,  then  the  new  pair  is  related  either  by 
case  a  or  b  (cf.  Figure  4)  according  to  whether  X  is  white  or 
black  respectively. 


Procedure  FINDNEIGHBOR  must  be  called  for  each  boundary 
segment  defined  by  a  (black, white)  node  pair.  Consider,  for 
example,  the  call  X* FIND  NEIGHBOR (P ,T, U) .  The  time  required 
is  equal  to  the  sum  of  the  path  lengths  from  nodes  P  and  X  to 
their  nearest  common  incestor.  Thus  in  the  worst  case  a  single 
call  will  take  hei jht  of  the  quadtree  time.  We  now  analyze  the 
average  case. 

In  a  complete  quadtree  of  height  n  there  are  2^n-2n  leaf 
nodes  which  have  east  neighbors.  Of  these  neighbor  pairs  2n 
have  their  nearest  common  ancestor  at  level  n,  2*2n  at  level  n-1, 

...,2  *2  “2  at  level  . . ind  2  at  level  1.  For 

example,  in  Figure  5  nodes  corresponding  to  squares  1-8  have 
eastern  neighbors  l*-8'  and  nearest  common  ancestors  at  level  3; 
nodes  for  squares  9-24  have  eastern  neighbors  9* -24'  and  nearest 
common  ancestor  at  level  2.  (These  same  frequencies  also  occur 
for  each  of  the  other  sides  north,  south  and  west.)  Therefore, 
if  we  assume  that  P  and  X  are  both  in  le\w*l  0  and  P  is  equally 
likely  to  occur  at  any  of  the  2*n  |>ossible  positions,  then  the 
average  time  required  by  FIND  NEIGHBOR  is 

(  ( 2n  •  2  n)  ♦  ( 2n* 1  •  2  (n- 1 )  )  ♦  •••  ♦  ( 22n~ 1  -  2  i  >  ♦•••♦<  22n_1  •  2)  )/( 220 

n  i 

=  I  2i/<21) 

1-1 

-  (2n<-1-(n*2))/2n‘1 

< 
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That  is,  P  and  X’s  nearest  ccuunon  ancestor  is  on  the  average 
at  level  2. 

When  P  and  X  are  both  at  level  k,  k>0,  a  similar  argument 
shows  that  on  the  average  only  a  constant  amount  of  time  is 
required  for  each  call  to  FIND  NEIGHBOR.  In  general,  if  P  is 
at  level  j  and  X  is  at  level  k,  then  the  average  tine  is  equal 
to  |  j-k  |  . 

The  average  total  time  required  by  all  calls  to  FIND  NEIGHBOR, 
and  consequently  NEXTLINK ,  depends  on  the  expected  disparity 
in  the  sizes  at  leaf  nodes  which  are  adjacent  to  leaf  nodes  of 
the  opposite  type.  If  all  leaf  nodes  are  at  level  0,  for  example, 
then  the  average  total  time  to  traverse  the  entire  boundary  is 
O (number  of  black  boundary  nodes).  Further  analysis  would  require 
knowledge  of  expected  region  shapes  and  positions  and  is  beyond 
the  scope  of  this  paper. 


7 .  Concluding  remarks 

An  algorithm  has  been  presented  for  converting  the  quadtree 
representation  of  a  simply-connected  region  into  the  chain  code 
description  of  the  region's  boundary.  The  boundary  traversal 
time  is  linear  in  the  number  of  boundary  nodes  in  the  juadtree 
when  certain  plausible  assumptions  are  made  about  the  form  of 
the  quadtree. 

In  the  case  where  a  region  may  have  holes,  we  may  extend 
the  algorithm  by  simply  adding  a  quadtree  traversal  procedure 
which  systematically  visits  all  black  nodes  upon  completion  of 
the  first  boundary  following  sequence .  If  that  scan  discovers  a 
black  boundary  node  having  a  boundary  edge  which  was  not  marked 
by  the  boundary  follower,  then  the  scan  is  temporarily  interrupted 
so  that  the  boundary  of  which  it  is  a  part  can  be  followed. 

There  are  a  number  of  other  problems  concerned  with  quadtrees 
which  are  of  interest  in  establishing  their  applicability  for 
efficiently  operating  on  areal  data.  The  algorithm  qiven  here 
assumes  a  single  connected  region;  an  efficient  algorithm  for 
determining  the  connectivity  of  the  black  nodes  would  be  of  con¬ 
siderable  value.  Another  important  extension  of  this  work  is 
to  consider  the  converse  operation  of  constructing  the  quadtree 
from  the  chain  code  of  a  region;  this  will  be  the  subject  of  a 


subsequent  paper. 


I 


Figure  2.  Labeling  of  a  block's  four  quadrants  and 
four  sides. 


Figure  3.  Possible  overlap  relationships  between  the 

(black .white)  adjacent  node  pair  (P.O) .  The 
heavy  arrow  indicates  the  boundary  segment 
just  output. 


g .  h .  i . 


Figure  4.  Possible  configurations  of  P,  0  and  their  neighbor 
blocks  in  determining  the  next  (black .white)  pair. 
Arrows  indicate  the  boundary  segments  associated 
with  the  old  and  new  pairs. 
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There  has  been  recent  interest  in  the  use  of  quadtrees  to 
represent  regions  in  an  image.  It  thus  becomes  desirable  to 
develop  efficient  methods  of  conversion  between  quadtrees  and 
other  types  of  region  representations.  This  paper  presents  an 
algorithm  for  converting  from  quadtrees  to  a  simple  class 
of  boundary  codes. 
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